AD-A157  912 


¥  933 


THEORY  OF  CONSOLIDATION  OF  SOFT  SEDIMENTS,  I 

(being  the  first  Interim  Report  to  U.S. A.R.D.S.G. 
under  Contract  No.  DAJA^^85-M-01 58  and. 
Requisition  R&DH^3 J-EN-01 ) 


-  U.K. 


Introduction 


The  objectives  of  the  research  undertaken  are  to  provide  a  theory  of 
consolidation  of  fine-grained  soils  of  sufficient  generality  to  enable 
predictions  to  be  made  in  engineering  situations  where  the  sediment  is  so  soft 
that  allowance  must  be  made  for  changes  in  material  properties  during  the 
progress  of  consolidation,  where  large  strains  and  displacements  must  occur, 
account  being  taken  of  pore  water  flow  in  two  or  three  dimensions. 

The  study  is  proceeding  in  two  Stages.  The  first  Stage  (being  undertaken 
largely  by  Prof.  R.E.  Gibson)  involves  a  critical  examination  of  presently 
available  theories  of  consolidation  which  are  modified  and  extended  to  meet  the 
objectives  of  the  research.  The  aim  in  Stage  I  is  to  develop  mathematical  / 

equations  governing  the  consolidation  process  based  on  physical  assumptions  ( 
which  accord  with  the  known  behaviour  of  soft  sediments.  \ 

In  Stage  II  (Principal  Investigators:  Prof.  R.L.  Schiffman  and  Prof.  R.E. 

Gibson)  numerical  procedures  and  associated  algorithms  will  be  considered  based  1 
on  the  analysis  which  would  permit  prediction  of  the  distribution  and  time 
variation  within  the  medium  of  quantities  of  engineering  interest,  namely:  the 
pore  water  pressure,  the  effective  stresses  and  the  void  ratio.  This  may  not  be 
straightforward  and  analytical  modifications  to  the  Stage  I  equations  may  be  / 
needed  to  secure  computable  results  which  can  be  generated  economically  and 
efficiently .  / 

This  first  Interim  Report  is  devoted  to  a  synoptic  account  of  the  progress  / 
achieved  in  Stage  I.-.. 

2.  Theories  of  Consolidation:  a  brief  Survey  / 

The  mechanism  of  the  consolidation  process  in  fine-grained  soils  and  its 
expression  in  analytical  form  were  given  first  by  Terzaghi  in  1923  ('). 
theory  was  restricted  to  pore  water  flow  in  one  (Cartesian)  dimension  but  as  was 
recently  pointed  out  did  take  some  account  of  large  strains  (2).  It  was  aimed 
at  providing  a  means  for  predicting  the  progress  of  settlement  in  uniformly 
loaded  clay  layers.  Since  that  time  the  theory  of  one-dimensional  consolidation 
has  been  greatly  extended  to  take  account  of  non-linear  behaviour,  including  the 
variation  during  consolidation  of  the  coefficients  of  permeability  and 
compressibility,  and  other  effects  associated  with  large  strain  and  deformation 
(see,  for  example,  ',  -,  *)  and  this  development  was  made  possible  by  the 


intrinsic  simplicity  of  one-dimensional  deformation. 
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The  problems  of  consolidation  which  face  the  civil  engineer  are  rarely 
one- dimensional  and  a  need  soon  developed  for  a  more  general  formulation  of  the 
theory  which  would,  for  example,  allow  the  settlement  of  structures  to  be 
predicted  or  sand  drain  installations  to  be  rationally  designed:  both  cases 
where  soil  deformation  and  pore  water  flow  occur  in  three  dimensions.  Severe 
difficulties  were  soon  encountered  in  this  quest  and  it  was  found  that  only  by 
introducing  radical  simplifying  assumptions  could  any  progress  be  made.  This 
phase  of  development  is  associated  with  Terzaghi  and  Rendulic  -  both  engineers  - 

and  with  the  mathematician  Maurice  Biot.  The  nature  of  their  contributions 

reflect  their  difference  in  outlook  and  objective.  Biot  developed  a  rigorous 
and  self-consistent  theory  (*)  based  on  the  assumption  that  the  soil  skeleton 
behaves  as  a  porous  perfectly  elastic  medium,  although  this  is  only  a  crude 
approximation  to  the  known  behaviour  of  real  soils.  Terzaghi,  on  the  other 
hand,  tried  to  avoid  commitment  to  particular  constitutive  relations  between 

components  of  strain  and  effective  stress  and  he  succeeded,  but  only  at  the  cost 

of  implying  that  during  consolidation  (under  constant  loading)  all  components  of 
total  stress  remain  constant  everywhere  (7,8).  In  reviewing  these  theories  in 
1943  (9)  Terzaghi  remarks  "...  the  existing  theoretical  methods  for  dealing  with 
two-  and  three-dimensional  problems  of  the  consolidation  of  clay  under  load  are 
not  yet  ready  for  practical  application".  Owing  to  the  complexity  of  the 
problem  progress  was  slow,  but  by  1969  the  consequences  of  these  two  differing 
approaches  had  been  examined  in  a  sufficient  number  of  particular  cases  for  some 
general  conclusions  to  be  drawn  by  Schiffman  and  co-workers  (l0). 

In  all  these  studies  the  assumption  was  invariably  made,  either  tacitly  or 
explicitly,  that  during  consolidation  the  effective  stress-strain  relations  are 
linear,  the  coefficient  of  permeability  of  the  soil  remains  constant  and  that 
the  strains  and  deformation  of  t*-e  soil  are  small.  However,  these  factors  can 
have  an  important  influence  on  the  progress  of  consolidation,  in  particular  upon 
the  magnitude  and  changing  pattern  of  pore  water  pressure  distribution.  More 
recently  engineers  have  sought  to  remove  those  restrictions  from  the  theory,  but 
owing  to  the  complexity  of  the  general  case  have  confined  their  attention  to 
problems  of  one-dimensional  compression  and  flow  (see,  for  example,  ‘•,,7,1B). 
Contemporaneously  Biot's  work  has  been  extended  to  large  strains  and 
deformations  (l*,20,21,22),  but  the  "soil"  considered  remains  highly  idealized*. 

3.  The  Present  Study 

The  typical  problem  with  which  we  are  concerned  here  is  the  consolidation  under 
its  own  weight  of  a  mound  of  soft  saturated  clay  which  has  been  placed  beneath 
water.  It  has  been  assumed  that  the  strains  and  deformations  that  may  develop 
during  consolidation  are  large:  that  is,  account  must  be  taken  in  the  theory  of 
the  variation  of  compressibility  and  permeability  during  consolidation. 
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Attention  has  for  the  present  been  restricted  to  plane  two-dimensional 
pore-water  flow  and  one-dimensional  soil  deformation.  The  details  of  the 
analyses  are  given  in  Appendices  A  (Eulerian  description)  and  B  (Lagrangian 
description).  The  most  promising  approach  seems  at  present  to  lie  with  a 
Lagrangian  description  allied  with  the  void  ratio  (e)  as  the  dependent  variable. 
The  governing  equation  (B19,  B26)  which  results  is  non-linear  and  must  be  solved 
numerically  using  an  iterative  procedure. 

A  numerical  study  of  specific  cases  may  reveal  that  strains  and  deformations  are 
not  unduly  large,  in  which  case  alternative  approximations  to  those  adopted  may 
prove  to  be  more  appropriate.  This  is  likely  to  be  helpful  if  the  more  general 
problem  of  successive  accretion  of  material  on  an  existing  mound  is  to  be 
treated  economically,  although  the  present  theory  can,  in  fact,  be  used  to  cover 
this  case. 

4.  Extension  to  the  Present  Study 

It  is  intended  that  the  concluding  parts  of  Stage  I  of  the  Study  should  be 
devoted  to  an  examination  of  the  feasibility  of  extending  the  work  recorded  in 
Appendix  B  to  cover  the  cases  of  axially  symmetric  and  three-dimensional 
consolidation  within  the  context  of  the  assumptions  adopted  therein. 
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APPENDIX  A 


In  this  Appendix  we  consider  the  case  of  pore-water  flow  and  soil  skeleton 
deformation  restricted  to  two  dimensions  (the  x,z-plane)  and  we  seek  to  derive 
the  equations  governing  the  motion  of  these  phases  by  adopting  the  so-called 
Eulerian  scheme  of  description.  We  therefore  consider  the  soil  grains  and  water 
which  at  any  time  (t)  reside  within  the  rectangular  element  of  sides  (fix,  6z) 
located  at  the  point  (x,z). 

We  denote  by  (wx,  wz)  the  velocity  components  of  the  solid  phase  and  by  (vx,  vz) 
those  of  the  pore  water.  These  components  will  be  functions  of  the  three 
independent  variables  (x,z,t). 

1 .  Equations  of  Continuity 

The  volume  porosity  (volume  of  void  space  Vv  per  unit  of  bulk  volume  V)  we 
denote  by  n,  so  that 

n  =  Vv/V, 

while  the  area  porosity  (area  of  void  space  Av  per  unit  of  bulk  area  A)  we 
denote  by  na.  If  the  number  of  grains  occupying  the  rectangle  is  in  some  sense 
large,  it  can  be  shown  that 


na  *  n. 

The  rate  of  increase  of  weight  of  solids  within  (fix,  6z)  must  equal  the  net  rate 
of  flow  of  weight  across  the  four  faces  of  the  element  and  this  leads  us  to  the 
equation 


[0-n)ps]  +  [wx(1-n)ps]  +  [wz(1-n)ps]  -  0  (A1) 

where  p3  is  the  unit  weight  of  the  solids. 

Similarly,  for  the  pore-water 

ft  tnpf^  +  tx  tvxnPf]  +  fi  [vZnPf]  -  0  (A2) 

where  pf  is  the  unit  weight  of  the  fluid.  These  are  the  equations  of 
continuity. 

Since,  in  the  range  of  stress  with  which  we  are  concerned,  the  unit  weights 
ps,  pf  can  be  regarded  as  constants,  it  follows  from  (A1)  and  (A2)  by  addition, 
that 


2. 


Flow  Rule 


We  shall  assume  that  the  pore-water  moves  through  the  soil  skeleton  in 
accordance  with  Darcy’s  law,  due  account  being  taken  of  the  fact  that  it  is  the 
relative  velocity  between  the  phases  which  induces  a  drag  on  the  soil  skeleton. 
Accordingly  we  write 


n(vx  -  wx)  -  -  & 

1  x  XJ  Pf  3x 

(AH) 

n(vz  -  wz)  -  - 

z  z;  Pf  o z 

(A5) 

where  (kx,  kz)  are  the  coefficients  of  permeability  appropriate  to  the  (x,z) 
directions  and  u  is  an  excess  pore  water  pressure  which  is  defined  in  terms  of 
the  pore-water  pressure  (p)  by 

u  =  p+  pfz  +  const. 

(A6) 

where  the  positive  direction  of  z  is  against  gravity. 

3.  The  Governing  Equations 

If  we  set  (A4),  (A5)  in  (A3)  we  find  that 

3  r.  3ui  3  r.  2u-i  r 3wx  ^  IWzI 

3x  kx  3jc  +  3z  kz  3z  '  pf[3/+3^] 

(A7) 

and  it  can  be  seen  that  if  the  soil  skeleton  is  stationary 
familiar  equation  governing  the  excess  pore-water  pressure 

(wx=0,  wz=0)  then  a 
is  obtained,  namely 

3  [i,  3ui  3  r.  3ut  n 

2^  Lkx  3^  3z  ^kz  8zJ  °* 

(A8) 

It  is  worth  noting  that  just  after  dumping  a  mound,  wx  «  0  and  wz  =  0  almost 
everywhere  and  so  (A8)  holds  at  this  instant;  this  equation  again  holds  when 

consolidation  is  complete. 
Now,  from  (A1 ) 


1 

(1-n) 


rlB 

l3t 


+ 


3n 

ciX 
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(A9) 


and  eliminating  the  rate  of  soil  skeleton  dilatation  between  (A7)  and  (A9)  we 
find 


□2u  ♦  ~  [in(l-n)]  - 


0 


(A10 


where  the  operator 


D  ,  3  3  3 

Dt  3t  +  W*  3x  +  Wz  3z 

denotes  a  differentiation  following  the  motion  and 


□  2 


1-  [hx  i_]  +  [iSz  i.1 

3x  lpf  3xJ  3z  Lpf  3zJ 


In  order  to  proceed  further  the  relations  connecting  strain  increments  with  the 
effective  stresses  and  the  effective  stress  increments  must  be  assumed.  These 
relations,  together  with  the  equations  of  equilibrium,  would  allow  a  complete 
theory  to  be  developed.  However,  in  view  of  the  wide  variety  of  soil  types  and 
their  complicated  response  to  load  a  completely  general  approach  is  at  present 
not  feasible. 


We  shall  therefore  proceed  in  the  spirit  of  Terzaghi  and  Rendulic's  early  work 
and  seek  plausible  simplifications  which  will  permit  solutions  to  be  obtained  to 
problems  which  will  be  sufficiently  exact  for  engineering  purposes. 

iJ.  One-dimensional  compression 

We  shall  assume  here  that  the  development  of  a  mound  takes  place  in  a  number  of 
stages  each  consisting  of  the  sudden  dumping  of  soil  followed  by  a  period  of 
consolidation  during  which  no  further  loading  occurs. 

During  the  consolidation  it  will  be  assumed  that: 

(a)  lateral  displacement  of  the  soil  can  be  ignored. 

The  mound  will  therefore  consolidate  with  the  development  of  vertical  strains 
only  but  the  pore  water  flow  we  shall  not  assume  to  be  constrained  in  any  way. 

The  porosity  (or  void  ratio)  at  any  point  in  the  clay  will  depend  on  the  initial 
porosity  and  state  of  effective  stress  there  and  also  upon  the  subsequent 
development  of  the  components  of  effective  stress  as  consolidation  proceeds. 

The  detailed  effective  stress-strain  behaviour  of  soils  is  very  complicated  and 
to  obtain  an  engineering  solution  to  the  problem  we  introduce  the  following 
simplifying  assumptions: 

(b)  The  void  ratio  (e)  of  a  soft  saturated  clay  depends  to  a  good 
approximation  only  upon  the  major  principal  effective  stress  o,'  (the 
so-called  "American  Hypothesis";  see,  for  example,  refs.  1 *",s). 

(c)  The  major  principal  stress  is  (almost  everywhere)  equal  to  the  vertical 
stress  o2z  (see  ref.1*)  in  a  mound-like  structure. 

It  follows  from  these  two  assumptions  that 

e  -  e(o,')  -  e(oZz’).  (All) 

For  simplicity  of  exposition  we  restrict  our  remarks  to  the  case  where 
kx«kz«k(const) ;  the  more  general  problem  can  if  required  readily  be  treated 
later  in  detail.  Equation  (A10)  can  be  written  in  the  form: 


Since 


where 


Pf  (1+e)  do,'  Dt 

.  D  3  3 

Mhere  Dt  *  St  *  “23T 


°zz  ~  P 


“v*’p  -  jjf-T? 


(A13) 


(A14) 


k(l+e)  da,' 
v  pf  de 

which  is,  apart  from  the  factor  (1+e)  in  place  of  (1+e0),  Terzaghi's  coefficient 
of  one- dimensional  consolidation. 

5.  One-dimensional  consolidation 


It  is  worth  noting  that  in  cases  of  strictly  one-dimensional  compression  and 
pore  water  flow,  the  above  assumptions  (a),  (b)  and  (c)  are  satisfied  and 
equation  (A14) 


3  2p  _  Dp  Da 


is  exact. 

This  equation  governing  the  pore  water  pressure  has  a  structure  very  similar  to 
that  encountered  in  Terzaghi's  theory  and  indeed  his  equation  can  be  recovered 
by  replacing  the  differential  operator  D/Dt  by  3/3t.  From  ( A1 3 )  this  can  be 
seen  to  be  tantamount  to  ignoring,  for  example,  wr3p/3z  compared  with  3p/3t. 

There  is  no  a  priori  reason  to  suppose  that  this  can  be  justified.  Some 
indication  of  the  error  involved  can  be  found  from  (A3)  which  reduces  to 

nv2f  +  (1-n)vzs  «  0  (A16 

when  there  exists  a  plane  on  which  the  condition  vzf  -  v2s  -  0  persists.  (This 
is  so  in  the  oedometer  test:  the  mid-plane  with  two-way  drainage,  or  the  base 
if  there  is  no  flow  from  it).  It  then  follows  from  (A16)  and  (A5)  that 


k  3u 
z  pf  3z 

and  so,  for  example: 


3 


9 


i 

What  then  is  the  equation  governing  the  pore-water  pressure  which  is  exact 
within  the  context  of  the  assumptions  made  in  the  theory,  when  the  Eulerian 
description  is  used? 

We  commence  from  the  form  of  (A15)  which  allows  for  the  variation  of 
[  permeability,  namely 


1 _  3  r  k  3u i 

m^  3z  pf  Jz' 


Dp  _  Do 
Dt  ~  Dt 


where,  from  (A6): 

Dp 
Dt 


Also 


9u  +  k_  3u  ,au  _  \ 

at  +  pf  3z  oz  '  PfJ 


Dq; 

Dt 


3o 

at 


k  8u  3c 


pf  az  az 

but  vertical  equilibrium  requires  that 


(A18) 


(A19) 


(A20) 


^l22  =  ”  npf  -  C 1  -n ) ps  (A21 ) 

Setting  (A21)  in  (A 20)  and  subtracting  the  resulting  equation  from  (A19),  it  is 
found  that  (A18)  becomes 


1  o  rk  oU-i  _  J< f 2  k(ps~pf )  3u  ^u  oo7.7.  (a??) 

mv  3z  Lpf  3zJ  pfl3z;  ~  (1+e)pf  3z  'at  Dt  * 

Using  (A6)  a  similar  equation  can  be  derived  for  the  pore-water  pressure  (p). 

This  equation  is  highly  non-linear  even  in  the  "thin"  layer  case  (ps=Pf)»  and 
its  use  in  problems  of  large  displacement  and  strain  is  limited  for  the 
following  reasons: 

(i)  The  parameters  k  and  mv  are  related  to  the  void  ratio  e,  but  this  in  turn 
is  not  connected  to  u  (or  p)  in  a  straightforward  way; 

(ii)  The  vertical  total  stress  czz  and  its  time  derivative  aozz/3t  for  a  fixed 
value  of  z  varies  with  time  in  a  way  which  is  not  given  as  part  of  the 
data  of  the  problem  and  must  be  discovered  as  part  of  the  solution; 

(iii)  The  geometry  of  the  mound  boundary  on  which  a  condition  of  the  type  u«0 
persists  is,  again,  not  known  ab  initio  and  must  be  found  during  the 
solution  process. 

However,  one  distinct  advantage  of  working  in  terms  of  u  is  that  the  initial 
conditions  can  be  specified  in  terms  of  the  initial  state  of  total  stress  in  the 
mound,  this  allows  a  more  complete  and  exact  description  than  can  be  achieved 
when  the  void  ratio  or  porosity  is  used  as  the  Independent  variable  (see 
Appendix  B) . 


A  rather  more  promising  governing  equation  emerges  if  the  porosity  (n)  is  taken 
as  the  dependent  variable.  Commencing  from  (A10)  we  can  write 


Now 


L.  iL  (k— ) 

Pf  3z  1  3z' 


1 

(l-n) 


3u  _  3p 
3z  ~  3z 


Pf 


Pf 


so  that  using  (A21 ) 


3u 

3z 


(1-n) (ps-pf ) 


Setting  this  expression  into  (A23)  we  find,  after  some  algebra,  that 


(A23) 


(A2iJ) 


(A25) 


(^  -  -)  I;  UO-nH  |2 


(A26) 


which  is  essentially  equation  (23)  of  (*).  It  is  not  open  to  objections  (i)  and 
(ii)  above.  The  difficulty  (iii)  remains  but  it  has  been  claimed  that  this  can 
be  overcome  by  using  suitable  numerical  techniques.* 


*Lee  and  Sills  (ref.5)  use  a  numerical  technique  originating  from  work  by  Crank 
and  Gupta  (,l»12)  which  effectively  updates  the  position  of  the  moving  boundary 
as  the  solution  proceeds.  They  take  the  case  of  a  "thin"  soil  layer  where  the 
approximation  ps  -  pf  can  be  justified  on  physical  grounds,  and  compare  their 
numerical  solution  with  that  given  by  Gibson,  England  and  Hussey  (Ref.^O.  The 
agreement  is  apparently  excellent,  but  is  not  wholly  convincing  as  the 
coefficient  cp  in  ref.1*  is  taken  incorrectly  as  cp  -  cv/(1+e)2  instead  of 
cp  -  cv(1+e8)2/(1+e)2. 


APPENDIX  B 


In  this  Appendix  we  consider  the  physical  problem  discussed  in  Appendix  A,  in 
particular  the  case  of  two-dimensional  pore-water  flow  and  one-dimensional 
compression  in  a  submarine  clay  mound.  However,  here  we  use  the  Lagrangian 
scheme  of  description  and  examine  the  motion  of  an  element  of  soil  which 
contains  the  same  solids  throughout  its  history. 


For  reasons  mentioned  in  Appendix  A  we  shall  not  seek  to  derive  an  equation 
governing  the  pore-water  pressure  (or  an  excess  pore-water  pressure),  but  rather 
work  in  terms  of  the  void  ratio  (e)  or  the  porosity  (n).  The  same  physical 
assumptions  used  there  will  again  be  adopted  here.  The  argument  which  ensues 
follows  closely  that  used  in  ref."  and  we  shall  not  repeat  the  details  here. 

1 .  Equations  of  Continuity 

The  deformation  of  the  soil  skeleton  is  described  by  following  the  motion  of  an 
element  of  initial  area  (5a  x  6b)  lying  at  t  =  0  at  (a,b)*  which  during  its 
subsequent  motion  always  contains  the  same  solids.  The  location  of  a  material 
point  initially  at  (a,b)  is  given  subsequently  by 


5  -  5(a,b,t)  (B1 ) 

T)  =  n(a,b,t)  (B2) 

where,  by  definition, 

a  =  5(a,b,o)  (b3) 

b  =  n(a.b.O).  (B4 ) 


Since  we  shall  assume  that  the  motion  of  the  soil  skeleton  is  confined  to  the  a 
direction  it  follows  that  (B2)  may  be  replaced  by 

n  =  b.  (B2)bis 


The  equation  of  continuity  of  solids  takes  the  simple  form 


where 


hQ  =  n(a,b,0) 
e0  =  e(a ,b,0) 


are,  respectively,  the  initial  (t  =  0)  porosity  and  void  ratio  at  (a,b). 

As  the  element  translates  and  deforms  during  its  subsequent  motion,  pore-water 
moves  into  and  out  from  its  four  faces  and  the  equation  of  continuity  of  this 
(incompressible)  phase  is  found  to  be 


la  tn(va-wa) ]  +  ffc(nvb|f)  +  =  °* 


Flow  Rule 


Darcy's  Law  takes  the  form* 


,  .  k  3u 

n(v  wa)  =  -  —  — 

Pf  oE, 


h(v b-wb)  =  - 


k  3u 
Pf  3n 


with  this  description,  but  since 

3u  _  3u  3 £ 

9a  3£  9a’ 

3u  3u 
3b  9n 

and  w^  =  0 

it  follows  that 


,  >  3E  k  3u 

n(va-wa)  7?  “  '  —  — 


*  We  have  taken  ka  -  k(,  >  k  here,  but  the  extension  to  the  anisotropic  case  is 
trivial . 


The  excess  pore-water  pressure  (u)  is  related  to  the  pore-water  pressure  (p)  by 
u  =  p  +  +  const.  (B14 

3.  The  Governing  Equations 
Setting  (B12)  and  (B1 3)  in  (B6)  we  find 


8  rk  3u/ 1  +e0 1 1  +  3  r  r  1  +e  i k  Su-i  _  3  r  e  i 
3a  pf  3a  1+e  JJ  3b  ll1+e0jPf  3bJ  "  St  M+e0J 


(B1 


where  (B5)  has  been  used.  Using  now  (B5)  and  (B1^),  equation  (B15)  can  be 
expressed  in  terms  of  the  pore-water  pressure: 


3  ^  k  3 p ^  1  +e0 
Vertical  e^uiffbri  ~ 


)  ^  k] 


+e 

um  requires 


3  r k  Sp^ 1 +e  1 1  1  3e 

thlf  ^  =  "H+eJ  3t 


(B1 


|f  +  [npf  ♦  ( 1  — n ) p3 ]  |f  =  0,  (B17 

where  o  is  the  vertical  total  stress,  and  this  may  be  used  in  conjunction  with 
p  =  o  -  o '  (B1 ' 

which  defines  the  vertical  effective  stress  (a’),  to  eliminate  p  from  (B16). 


The  resulting  equation  governing  the  void  ratio  (e)  is 


3  r_Cp  3e-i  3  (1+e)2  3ei 

3aM+e0  3aJ  3bLL*  (1+e0) J  3bJ 


3  rk  rl+e  3 o -[ 
3blpfl1+e0 j3bJ 


1  3e 
( 1 +e0 )  3t 

in  the  general  case,  which  reduces  when  e0  is  a  constant  to 


fh[-(1+e)ff] 

ob  pf  3b 


at 

(BIS 


where  (as  in  ref.4)  the  coefficient  of  finite  consolidation 


The  total  stress  term 


4. 


The  last  term  on  the  left-hand  side  of  (B19)  is  not  known  a  priori  and  must  be 
found  and  updated  during  the  solution.  In  this  Section  we  show  how  the  term 
3o/3b  can  be  found  at  any  time  from  the  current  void  ratio  distribution 
e(a,b,t) . 

Consider  a  vertical  cylinder  of  material,  of  unit  cross-sectional  area, 
extending  from  the  base  of  the  mound  to  the  water  level  above  the  mound  and 
denote  this  height  by  H.  It  is  evident  that  the  weight  of  material  (solids  and 
water)  within  this  cylinder  would  remain  constant  if  no  water  flowed  across  the 
surfaces  (i.e.  strict  one-dimensional  consolidation)  or  was  added  (H  -=  constant). 
Under  these  conditions  3a/3b  on  the  base  depends  on  b  alone.  However,  on  the 
surface  of  the  mound  a  =  a0(b)  settlement  increases  the  water  pressure  and 
hence  the  total  vertical  stress  there.  Thus  at  the  top  of  the  clay  3o/3b 
depends  both  on  b  and  t_.  At  any  other  height  a,  or  when  the  conditions  of 
simple  vertical  flow  of  pore-water  no  longer  obtain  or  H  varies  with  time,  the 
term  3o/3b  can  be  expected  to  depend  upon  (a,b,t).  We  shall  now  determine  this 
relation. 

Integrating  (B17)  with  respect  to  a,  and  using  (B5),  it  is  found  that 


a  "  pf  r°(b)TTi7da  +  ps  r(b>  TtJ;  +  p^H(t)  *  C(a.(b),b.t)]  (B 22) 

where  the  last  term  represents  the  water  pressure  on  the  surface  a  -  a0(b)  of 
the  mound,  which  equals  the  total  stress  o  in  (B22)  when  a  »  a0(b). 

From  (B5): 


5(a,b,t)  =  J  j 


a  1+e(a,b,t) 
+e0 (a,b) 


da 


(B23) 


since  5(0, b,t)  *  0.  Substituting  (B23)  into  (B 22)  and  re-arranging,  it  is  found 
that* 

0  -  pfH  -  pf£a  (^)da  +  (Ps'Pf)^0  -f— 

It  follows  from  (B24)  that 


(B24) 


is  -  -  >ff  a. .  r- 


0 

(B25) 


*  The  first  term  represents  the  weight  of  water  if  it  completely  filled  the 
cylinder  from  the  base  to  the  water  surface.  The  second  term  is  the  weight  of 
water  if  it  filled  the  cylinder  from  the  base  to  the  height  £(a,b,t)  at  which 
0  is  acting.  The  last  term  is  the  buoyant  weight  of  solids  above  the  plane  on 
which  0  acts. 


const,  to 


which  reduces  in  the  case  e0  *= 


So  _  (ps~pf )  da 
8b  “  1+e0  db 


da , 


( 


the  first  term  arising  from  the  initial  non-uniform  surface  profile  of  the 
mound,  the  second  being  a  time- varying  correction  to  this  resulting  from  the 
settlement. 


5.  Some  typical  boundary  conditions 


If  a  mass  of  saturated  clay  of  more-or-less  uniform  void  ratio  e0  is  deposited 
in  a  short  period  of  time  beneath  water  and  rests  on  (say)  an  impervious  base, 
the  void  ratio  distribution  in  the  mound  at  subsequent  times  will,  on  the  basis 
of  the  assumptions  that  have  been  made  herein,  be  governed  by  (B19)  with  (B26) 


and  subject  to: 

(a)  the  initial  condition 

e(a,b,c)  =  e0(constant);  (B27 

(b)  the  mound  surface  boundary  condition 

e(a0 (b) ,b,t)  -  e0 (constant) ;  (B28 

(c)  the  impervious  base  condition 

wa(o,b,t)  =  0  (B29 

va(o,b,t)  -  0  (B30 


which  result  in  the  following  condition  on  the  void  ratio: 

9e  +  (pR-pf )  de 
8a  (1+e0)  do' 
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(B31 
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